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ABSTRACT 


A theoretical study and computer analysis of the cathode 
spot of a unipolar arc was conducted. The underlying theories 
of Plasma Physics, Space Charge Effects, and Electron Emission 
needed for an understanding of cathode phenomena are presented. 
Two models of unipolar arcing are reviewed and an analysis of 
the cathode spot is begun. A stationary model of the cathode 
spot is formulated with a system of equations that is not 
closed. A method of solution, using the Steenbeck Minimum 
Principle, is utilized and a computer program developed to 
determine the cathode spot parameters. Results of the analysis 
are presented for arc currents of 100, 150, 175, 200, 300, and 
400 Amperes. It is found that for copper stable arcs occur 


only for arc currents greater than 200 Amperes. 
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I. INTRODUCTION 


The study of cathode phenomena in gaseous discharges and 
electrical arcs has progressed over the last hundred years. In 
the last twenty years, this research has been revived as the 
interest in high power devices for pulsed power applications, 
Switching devices, and nuclear fusion power generation has 
increased. The discovery in 1958 of a new form of electrical 
arc, one that burns between a metal plate and a dense plasma, 
called a unipolar arc has further increased interest in the 
cathode phenomena. For both the bipolar and unipolar arcs, the 
current transfer occurs in localized regions between a cathode 
and plasma which is referred to as the cathode spot. It has 
become increasingly apparent that the phenomena occurring in 
thew cavmoce Spot Of a unipolar arc and that occurring in a 
bipolar arc are similar. Therefore, much of the theoretical work 
done on the cathode spot of a bipolar arc can be appled 
towards the understanding and formulation of a consistent 
theory of the unipolar arc. 

This study will begin with a short review of some basic 
theories of plasma physics, space charge effects, and electron 
emission mechanisms. Two useful unipolar arcing models are then 
reviewed. The stationary cathode spot will then be studied and 


a set of equations presented that can be used to determine it’s 


parameters. It will be found that the system of equations is 
not closed, hovevem a method for the solution of the system 
can be formulated by using the Steenbeck Minimum Principle. 
Finally, results of calculations based on these equations is 


presented. 


10 


If. BACKGROUND THERORY 


Ae PLASMA PHYSICS 
L Plasma Characteristics 

A plasma is formed when atoms are ionized into positive 
ions and negative electrons such that they form ae gas. 
However, not all such gases are considered to be plasmas. A 
more complete definition of a plasma has been given by Chen 
fRef...1); 

"A plasma is a quasineutral gas of charged and neutral 

particles which exhibits collective behavior" 

“Collective behavior’ means that the behavior of a small region 
of the plasma is determined not only by conditions in the plasma 
that immediately surrounds it, but also by conditions in the 
plasma far away from it. 

A plasma can be characterized by two quantities, n, the 
plasma density and T, the plasma temperature. Since the plasma 
must be quasineutral the density of ions, ni, and the density of 
electrons, me, must be approximately equal Therefore, the 
plasma density is defined as 

nny, = Ry. ( eqneezy! ) 
The plasma temperature is generally given in units of energy 
which corresponds to kT, where k is Boltzgmann’s constant and 7 


is the plasma temperature in Kelvin. For a gas in thermal 


te 


equilibrium, which follows a Maxwellian distribution, the average 
energy of the gas particles is related to the plasma 


temperature by 


(eqn. 2:2 
E a= SkT 


where the plasma is assumed to be three dimensional 
2. Debye Shielding 

When a plasma is subjected to an external electric field 
the electrons and ions will drift in opposite directions. As 
they drift apart, they produce a charge separation and an 
internal electric field, which opposes the external field. The 
width of the region, over which the charge separation must occur 
to balance the external electric field, is proportional to the 
thermal energy of the plasma particles. If the dimension of the 
bulk plasma is greater than the dimension of the region over 
which this charge separation occurs, the interior of the plasma 
will be shielded from the external field. 

The width of the charge separation can be determined, 
approximately, by considering a plasma where the ions are fixed 
in space, over the time frame of interest, and the electrons 
obey a Maxwellian distribution such that the Boltzman relation 
applies. Using Poisson’s Equation in one dimension gives 

d*Vv e C eqn. 2.3 ) 


2 
dx Eo 
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where V is the electric potential and e is the electric charge. 
Since the ions are considered fixed, the density of ions will be 
constant and equal to the plasma density no. The Boltzman 
relation then gives the electron distribution in the region of 


interest as 





eV ( eqn. 2.4 ) 
Ne=NoeXP| 7 |. 
8 


Substituting equation 2.4 into equation 2.3 gives 


a4 an eV 1 ( eqnee@es ) 
a a exp kT. 


To solve equation 2.5 the exponential can be expanded in a 








Taylor series, and neglecting terms of second order and higher 
gives 


aay ne2 @eqn 276 ) 
dx? EgkT, 





The solution of equation 2.6, which is a homogenous second order 


linear ordinary differential equation is given by 


(eqns 7 -) 
VY = v oexp( -*] 


D 


ak le ( eqn. 2.8 ) 
Ap=4/——- 
re 


This quantity is known as the Debye Length and is a measure of 


where 


the thickness of the shield that screens the bulk plasma from 
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the affect of external fields. This screening applies as long as 
the ae enone of the bulk plasma is much greater than the 
Debye Length. 
3. Sheaths 

When a plasma comes in contact with a wall the electrons 
and ions that hit the wall will recombine. But since the 
electrons are moving at a higher velocity than the ions, the 
electrons will be lost faster than the ions resulting in a net 
positive charge at the plasma-wall boundary. The wall will then 
be at a lower potential than the bulk plasma and an electric 
field will exist. Figure 2.1 illustrates the potential variation in 
a plasma which is in contact with a wall Due to Debye Shielding 
a layer of charge separation will exist, called the sheath, which 
will isolate the bulk plasma from this potential difference. The 
effect of the sheath is to accelerate the ions and to 
decelerate the electrons that enter the region until the flux of 
ions is balanced by the flux of electrons. Therefore, only those 
electrons in the high energy tail of the velocity distribution 
will be energetic enough to cross the sheath and all others will 
be repelled back into the plasma. 

The variation of the potential in the plasma sheath will 
now be considered. In order to simplify the problem the 
following assumptions will be made: 


* The ions enter the sheath region with a drift velocity 
Mon 


* The ions have Ti=0 so that all ions have a velocity uo 
at the plasma-sheath boundary, 
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* The sheath region is collisionless and in steady state, 
* The potential decreases monotonically with x, 


* The electrons follow the Boltzman relation (equation 
2.4). 





Figure 2.1 Plasma-Wall Interaction 


Applying conservation of energy to the ions in the sheath 
region gives 


=Mu?==Mus-eV (x). ( eqn. 2.9 ) 


Solving equation 2.9 for the ion velocity gives 


us ug- ( eqn. 2.10 ) 





The ion density within the sheath is determined from the ion 
equation of continuity which is given by 


Mise IL XL 5 Uo - (ean, 2.41.) 


i 


Using equation 2.10 for the ion velocity, equation 2.11 can be 


solved for the ion density which gives 


2eVv 
n(x) = Ro 1 = 2 
uoM 
Substituting equation 2.4 and 2.12 into equation 2.3 gives 


ot 
d*V eo (ev). ;_2eVv \ ° 
dx? Eo © kT, uzM 


Equation 2.13 is the non-linear planar sheath equation. By 


( eqn. 2125) 


Ni = 





{ eqn. 27 hom 














multiplying each side by the first derivative of the potential 


the equation can be integrated once to give 
V 2S 
. 
2NokT, eV Mus 2ev \* 
———*( | exp -1|+ 1-—— ] -1 ] ). 
Eo kT, kT, Mus 


( eqn. 2.14 ) 











Further solution of equation 2.14 would require numerical 
methods. If the electric field inside the plasma is zero, then 
the first derivative of the potential inside the plasma must 
also be zero. The left hand side of equation 2.14 is therefore 


greater than zero, and the following inequality results 


re > ( eqn. 2.15 ) 
u 
exp{ 22-1 2 (1-284) - 1 > 0O. 











kT, kT, Mué 
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This inequality can be greatly simplified by expanding the left 
hand terms in Taylor series expansions and neglecting terms of 
third order and higher. The resulting inequality is known as the 


Bohm sheath criterion and is given by 


ele.. ( eqn. 2.16 ) 
Uy ? M : 


This requires that the ions must enter the sheath region with a 





velocity greater than the ion acoustic velocity of the plasma. 
In order for this to occur, there must be a finite electric field 
within the plasma. Therefore, the assumption that the potential 
and the electric field are zero at the plasma-sheath boundary is 


only an approximation. 


B. SPACK CHARGE LIMITED CURRENT FLOW 
L Child-Langmuir Law 

The current that can flow between the cathode and anode 
of a vacuum diode is limited by the space charge of the 
electrons that exist between them. The electrons distort the 
external field and thereby reduce and can even reverse the field 
at the cathode surface. 

As an example, consider a diode consisting of infinite flat 
plates separated by a distance da Assume that there is an 
unlimited supply of electrons available from the cathode, and 
that the initial velocity of the electrons after emission is 


negligible compared with the velocity that they will gain while 
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crossing the electrode gap. The kinetic energy of the electrons 
crossing the gap and the current density carried by the 


electrons 1S given by 


=mu?= eV ( eqn. 2.17 ) 


and 
j=enu. (ean. 2.tem) 
Using Poisson’s equation, in one dimension, the current-voltage 


relationship can be derived as follows: 


aay ee) _ en(x) 


3 - (eqn. 2.19 ) 


The electron velocity is determined from equation 2.17 to be 


2e 


u( x)= ar ( eqn. 2.20 ) 


The electron density is determined using equation 2.18 and 2.20, 


and after substituting into equation 2.19 gives 


ayo 3 
cag ee ae 
dx? Eo 2e 





("@aqn 2.2 i 


In order to solve equation 2.21 both sides are multiplied by the 
first derivative of the potential and reduced to the following 


form 


d [dV \*inj [2a 
—| — |] =—,/ —V *—. , oe 
2() d eed e 
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Equation 2.22 is then integrated over the gap using the 
boundary conditions that the electric field and the potential are 


zero at the cathode. The result of this integration gives 


av _ [2i [am ae & 
dx E> a poem.) 


Integrating once again over the gap, applying the boundary 
conditions, and rearranging gives the desired current-voltage 
relationship of 


 4y¥2€, eV? 
J 9 mx? 


( eqn. 2.24 ) 

Equation 2.24 is known as the Child-Langmuir Law for 
space charge limited current flow [Ref. 2]. Although it has been 
derived here using simple assumptions, the proportionality of the 
current density to the 3/2 power of the potential difference 
remains for more difficult geometries under non-relativistic 
conditions. It should also be pointed out that the Child-Lang- 
muir Law applies to both positive and negative charge carriers. 
Figure 2.2 is a plot of electron current density versus 
potential for several gap widths. 

2. Space Charge Limited Current in a Sheath 

In deriving the planar sheath equation it was assumed 

that the electron density followed the Boltzman relation 


equation 2.4. Therefore, in a region close to the wall the 


electron density will be much less than the ion density 


i 
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and can be neglected. The first derivative of the potential at 
the point in the sheath where the electron density can be 
neglected, is also assumed to be negligible compared with the 


slope of the potential near the wall Equation 2.13 then reduces 














TO. 
eet | ( eqn. 2.25 ) 
ace Eo uZM eet Ss 
Making the change of variable that V = -V and noting that 
2eV 
>» | eqn. 2.26 
u2M ( eq ) 
equation 2.25 reduces to 
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The ion current density is constant across the sheath and is 
given by 

J/c=ONoto. Geqn. 2.25.) 
After substituting equation 2.28 into equation 2.27, the 
resulting expression is the same as equation 2.21. Therefore, 
the solution of equation 2.27 for the ion current density 
results in the Child-Langmuir Law for space charge Ilmited 


current flow, equation 2.24. 


Za 


3. MacKeown Equation 

In = arc, the current is carried primarily by the 
electrons emitted by the cathode and by a smaller amount of 
positive ions flowing from the plasma to the cathode. If the 
Space charge, which causes the cathode drop in the plasma 
sheath, is high enough it may exert a strong electric field at 
the surface of the cathode. This can then result in a large 
field emission current in addition to any thermionic current. 

The electric field at the cathode surface can be 
determined using Poissons Equation The following derivation 
measures distance and potential difference from inside the 
plasma where the cathode drop begins. The electron and ion 
current densities are given by 


Je=en,v C eqn. 2.29 ) 


and 
J, =enu. ( eqn. 2.c0)) 
Substituting equation 2.29 and 2.30 into equation 2.3 gives 


at --i(4-4] (eqn. 2sa 1m) 


dx? Eo U v 





where V is the potential difference, ji, the ion current density, 
je, the electron current density, u, the ion velocity, and v, the 
electron velocity. 

The velocity of the ions and electrons are determined 


using the conservation of energy. The initial energies are 


ee 


neglected, based on the assumption that the energy gained by 
the ions and electrons in crossing the sheath are much greater 


than the initial energies. For the ions the energy balance is 


] ( eqn. 2.32 ) 
—~Mu*=eV 
2 
and for the electrons the energy balance gives 
] (eqns 2-33) 


a 2_ _ 
Bee e(V. V } 


where M is the mass of the positive ion, m, the mass of the 
electron and Vc is the potential of the cathode. Substituting 


equation 2.32 and 2.33 into equation 2.31 gives 


< ( eqn. 2.34 ) 


1 
dave 1) (6) -1 m : 
dx? es} -‘\2ev % 2e(V.-V) 


Multiplying both sides by the first derivative of the potential 





difference, equation 2.30 can be integrated once and gives 


mts (ey Ae ae Va) ie, mV ,\? 
» Es ioe) 2e Je\ “De 


(eqn. 2.33 ) 





where F is the electric field. Evaluating for V =v gives the 


electric field at the cathode surface as 


1 ( eqn. 2.36 ) 


| j 1-2( 2) 
Te | j\M 





ae 


Equation 2.36 is known as Mackeown’s Equation [Ref. 3] and 
can be used to determine the electric field existing at the 
cathode provided the values of the cathode drop, the ion 
current density, and the electron current density are Iknown. 
Figure 2.3 is a plot of electric field at the cathode surface 
versus ion current density for several ratios of electron to 


ion current densities. 


ce ELECTRON EMISSION FROM METALS 

The theory of electron emission from metals is very well 
understood. The basic processes have been explained on the 
basis of a free electron gas in the metal, that is described by 
Fermi-Dirac statistics, and a potential barrier at the metal 
surface. Figure 2.4 shows a typical potential energy diagram 
for an electron as a function of distance from the metal-vacuum 
interface. For thermionic emission, that is emission at high 
temperature and low electric field, the electrons are emitted 
predominately by passing over the potential barrier at the metal 
surface. For field emission, that is emission at low temperature 
and high electric field, the electrons are emitted predominately 


by tunneling through the potential barrier. 
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Figure 2.3 MacKeown’s Equation for Electric Field at the 
Cathode Surface 
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Figure 2.4 Potential Energy Diagram for a Metal-Vacuum 
Interface 


In the past, thermionic and field emission have been treated 
separately, which has resulted in the well kmown Richardson 
equation for thermionic emission, and the Fowler-Nordheim 
equation for field emission The expressions that will be used 
here, however, are based on an approach that treats both 
thermionic and field emission from a unified point of view 
(Ref. 4}. 

Hartree units are used in the following sections because 
they greatly simplify the expressions. Thus the following 


quantities are redefined as: 


6:2 eT Ole 
j mi (eqn. 2.37) 
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where the starred quantities on the left are in SI units and 
their respective counterparts on the right are in Hartree units. 
The following sections will give expressions that can be used 
to determine the current density that can be expected to be 
emitted from a metal given the surface temperature, the electric 
field at the surface, and the work function for the metal under 
consideration. In addition, a set of applicability conditions are 
given that can be used to determine the validity of the 
particular expression for a given situation. Figure ZA 
illustrates the regions of applicability of the emission 
equations on a plot of temperature versus electric field. 
Appendix A is a Fortran program that calculates the current 
density that can be expected to be emitted. Appendix B and C 
are Fortran subroutines that calculate the current density and 
applicability of the expressions to be used respectively. 
L Thermionic-Field Emission 
The current density that can be expected from thermionic 


emission, corrected for the Schottky effect, is given by 


a 
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Figure 2.5 Regions of Applicability for Thermionic, Field, and 
Intermediate Electron Emission Expressions 


For equation 2.41 to be applicable the following conditions must 
be satisfied: 


sae} 1 a “( _ : ( eqn. 2.43 ) 
in{ d ata) = 


in{ A=S)) l — G@eaqn, 2.44 ) 
‘d (l1-d} ; 


Equation 2.41 reduces to the Richardson equation in the limit as 
the electric field is reduced toward zero. Appendix D is a 
lis ting of a Fortran subroutine that calculates the 
Thermionic-Field emission current density based on equation 2.41. 

2. Ejeld-Thermionic Kmission 
The current density that can be expected from field 


emission, corrected for temperature, is given by 


F? ( mtckT | 4y2¢7v 
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Equations 2.46, 2.48, and 2.49 are the Fowler-Nordheim Elliptic 
functions. The functions K[k] and E[{k] are complete elliptic 
integrals of the first and second kind respectively and are 


defined as 


Ktk}= f (1 -k?sin2@) *a@, (eqn. 2.51 ) 


E{k)= f (1 -k?sin2e)'a@. (eqn. 2.52 ) 


For equation 2.45 to be applicable the following conditions must 


be satisfied: 





- pf kT 
- F?> — : ; 
¢-F ae ate ( eqn. 2.0cm) 
l-ckT >V2fkT, ( eqn. 2.54 ) 
(29 
= 2198 
f DF(¢2-F)”? ( eqn ) 


Equation 2.45 reduces to the Fowler-Nordheim equation in the 
limit as the temperature approaches zero. Appendix E is a 
Se ine of a Fortran subroutine that calculates the 
Field-Thermionic current density based on equation 2.45. This 
subroutine makes calls to the subroutines listed in Appendix F, 
G, and H, which are used to calculate the Complete Elliptic 
Integral of the First Kind, Complete Elliptic Integral of the 
Second _ Kind, and the Fowler-Nordheim Elliptic functions 


respectively. 


30 


Seeicd liate Field-T) Po Wane 33 
When equations 2.41 and 2.45 are not applicable, the 

current density from field-thermionic emission may be determined 

in the intermediate region as 

Perera $ F26 

f= ——}| — | exp| - — + ——— Bee. 

J -() o( or aS C eq 


where 


9-5-5. ( eqn. 2.57 ) 


The argument of the functions t, equation 2.46, and v, equations 


2.48 and 2.49, as used in equations 2.56 and 2.57, is given by 


ed 


Bi 
Gi 


Eta can be determined by solving the folloWing equation by 
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For equation 2.56 to be applicable the following conditions must 


be satisfied: 
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Appendix I is a listing of a Fortran subroutine that 
calculates ne intermediate Field-Thermionic emission current 
density based on equation 2.56. 

4. Explosive Electron Emission 

The emission mechanisms discussed so far have wide 
applicability to many physical devices. However, with the advent 
of high current cathodes the thermionic and field emission 
theories were unable to explain the high current densities 
observed. One explanation for this high current density is that 
the surface of the cathode is not smooth, but actually has a 
surface covered with micropoints, or whiskers. The significance 
of the whisker is that the electric field, at the surface of the 
whisker is enhanced by a factor of 10 - 1000. Figure 2.6 
illustrates the electric field enhancement that can occur at the 
tip of a whisker. This results in a higher field emission 
current density from the whisker, and ultimately leads to 
explosive destruction of the emission site. This destruction 
produces a local burst of plasma that expands rapidly away from 
the cathode surface and has therefore come to be known as a 
cathode flare or plasma jet. The electrons, that reach the 
anode, are then emitted from the surface of the plasma jet. 
This whole process is commonly referred to as explosive electron 


emission (Refs. 5 and 6]. 
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Figure 2.6 Electric Field Enhancement at a Whisker 


Fursei and Zhukov [{Refs. 7, 8, and 9] have experimentally 
characterized the explosive emission process and have confirmed 
the existence of micropoints and that the current density 
follows the Fowler-Nordheim field emission theory with the 


inclusion of a field enhancement factor. 
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Itt. UNITPOLAR ARCING 


A. Robson-Thonemann Model 

In 1958, Robson and Thoneman observed an electrical arc 
burning between a dense hot plasma and a metal wall {Ref. 10]. 
What made this unusual was that there was no anode present to 
collect the electrons emitted from the conducting wall 
Apparently, the wall was acting as both a cathode and anode, and 
therefore, they named what they observed a “unipolar arc.’ 

Robson and Thoneman proposed the following model to explain 
the occurrence of unipolar arcing. When a hot dense plasma 
comes in contact with a conducting wall a sheath will form [If 
the plasma temperature is high, the sheath or floating potential 
can exceed that necessary to initiate an electric arc A 
cathode spot, which is a local center of emission, could then 
form on the metal surface. The large electron current flow into 
the plasma from the arc would result in lowering the plasma 
potential to the cathode fall potential of the arc. This allows 
more electrons from the high energy tail of the velocity 
distribution to cross the sheath and recombine at the wall and 
thus closing the current loop. Figure 3.1 illustrates’ the 
situation in the sheath prior to initiation of an arc, and Figure 


3.2 shows the return current flow after arc initiation 
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Figure 3.1 Sheath Region Prior to Unipolar Arc Initiation 
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Figure 3.2 Sheath Region After Unipolar Arc Initiation 


The resulting change in current flow outside of the arc is 


then given by 


od ie, Ae eV. eV, 
Meee) lo \ RTA) NRPS 











SD 


(eqn. 3:1.) 


Where I is the current, A is the area of the metal plate, and Ve 
is the pathos drop of the arc. But since the net current flow 
between the plasma and the wall must be zero, equation 3.1 also 
gives the current flow inside the arc. 

The ion current density in the sheath can be approximated 
as 
kT. ( eqn. 3.2 ) 
M 


where the Bohm sheath criterion is used to approximate the ion 





J, 7 eng 


velocity. The electron current density in the sheath can be 
determined by noting that the net flux of electrons, * from the 
plasma into the sheath is given by 

( eqn. 3ocm 


Lia 
b=—n,v. 
qed 


Using the Boltzman relation for the electron density and the 
expression for the average speed of particles following a 


Maxwellian distribution gives 


] eV, SkT, (eqn. 3.4 ) 
Ja= [eR , exp ae 
4 kT, wm 


The net current density across the sheath must be zero. 





Therefore, setting equation 3.2 equal to equation 3.4, and 


solving for the sheath potential gives 


kT, ( M ( eqn. 3.5 ) 
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Robson and Thoneman point out that there is a minimum 
Current necessary to maintain an arc burning. Therefore, 
according to equation 3.1, if the plate area is small enough a 
unipolar arc can not occur. For example, using copper with a 
plasma temperature of 4 eV, a cathode fall potential of 15 V, a 
plasma density of 1013 cm-3, and an are current of 100 amperes 
would require a plate area of 115 cm2. For a circular area, the 


radius required is 6.05 cm. 


B. Schwirzke-Taylor Model 

Schwirzke and Taylor have expanded on the basic concepts of 
the Robson-Thoneman model [{Ref. 11). They propose that for an 
arc to occur, the ion density above the cathode spot must 
increase to allow for the arc current to flow. This is contrary 
to the Robson-Thoneman model which assumes a constant plasma 
density. Once a cathode spot forms the desorption of neutral 
atoms and evaporation of the metal atoms from the spot is 
greatly increased. If only a small fraction of these neutral 
atoms become ionized in the sheath the plasma density will 
increase above the cathode spot. Thus, according to equation 
2.8 the Debye length will decrease resulting in a decrease of 
the sheath width above the cathode spot. The electric field 
above the spot increases, where the field can be approximated 


as 





V, Cean. 3.0") 
| mun ,kT in Bia 
D 


af 


Schwirzke (Refs. 12 and 13] elaborates further on the 
electric fields set up inside the high density Plasma which forms 
above the spot. The fluid equation of motion of the electrons 


in the plasma is given by 


dv 


ns rs 2 oe en,. ( eqn. G7 ) 
nsm—"=-en,(F+u,xB)-VP,+ ; 


J 





The time derivative can be neglected over the time span of an 
arc, and further assuming no magnetic field, equation 3.7 can be 


solved for the electric field, which gives 


j l ( eqn. 3.8 ) 





The first term on the right hand side of equation 3.8 can be 
neglected in a weakly ionized plasma. This implies that the mean 
free path length for collision between the electrons and the 
neutral atoms is much greater than the dimensions of the region 
of increased density above the cathode spot. Thus, the 
increased plasma pressure above the spot creates a radial 


electric field given by 


( eqn. 3.9 ) 


= l kT, 
en, en, 








Vine. 


The affect of this radial electric field is to lower the plasma 
potential in a ring around the cathode spot. The change in the 
potential can be determined by integrating equation 3.9 from the 
center of the arc, where the plasma density is highest, radially 


outward until the plasma density equals the background density. 


oO 


The change in potential is then given by 


nice. ie Geqnms.10 ) 
AV(r})=- In 
e Roo 








In this ring of reduced potential, electrons in the tail of the 
Maxwellian distribution can cross the sheath more easily and 
thus a higher return current flow from the plasma to the wall 
results. This current flow closes the current loop of the 
unipolar are. The concepts of the Schwirzke-Taylor model are 


illustrated in Figure 3.3. 
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Figure 3.3. Schwirzke-Taylor Unipolar Are Model 
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If equation 3.10 is set equal to the floating sheath 
potential equation 3.5, the ratio of the electron density in the 
plasma above the spot to the background plasma density which 


reduces the sheath potential to zero can be determined as 


Ng - M ( eqn. 3.11 ) 
Ro. YeQuann | 


Equation 3.3 and equation 3.11 can be used to calculate the 





maximum area necessary to support a given unipolar arc return 
current. By setting the current equal to the current derived 
for the Robson-Thoneman model, equation 3.1, the ratio of the 
return current areas required by each model can be compared. 


The result of the comparison gives 


Ag-r _ M eV. 2nm ee M (eqn. 3.12 ) 
As., V2um\~~P\ kT.) VOM Vi 2num 


where AR-T and AS-T are the areas required by the 








Robson-Thoneman model and Schwirzke-Taylor model respectively. 
As an example, for copper with a cathode fall potential of 15 
volts and a plasma temperature of 4 eV the area ratio is 8.43 x 
103. Therefore, for a 100 ampere arc the area required by the 
Schwirzke-Taylor model is only 137 x 10-2 cm2, or for a circular 


area the radius required is only 6.59 x 10-2 cm. 
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IV. CATHODE SPOTS 


A. TYPES OF CATHODE SPOTS 

A cathode spot is a small, luminous region on the cathode, 
where a localized transfer of current occurs between the 
cathode and a plasma in an arc discharge. The type of arc 
discharge can be either bipolar or unipolar, however in all cases 
a cathode spot is present. In addition, the spots are 
associated with the erosion of the cathode and therefore, are 
also localized centers for the transfer of material from the 
cathode into the plasma. This has been evidenced by the 
formation of craters on the surface of the cathode. 

It has been possible over many observations to distinguish 
between two different types of cathode spots [Ref. 14]. The 
first type, or type I spots, are characterized as rapidly moving 
spots with speeds of approximately 104 cm/s across the surface 
of the cathode and with relatively low erosion from the surface. 
The type II spots, on the other hand, have speeds of 
approximately 102 cm/s and have much greater erosion from the 
surface. The type Il spot also has a tendency to group together 
and form clusters. The mechanism of erosion for the type Il 
spot is believed to be due to non-thermal means, such as 
explosions of whiskers or microinhomogenieties on the cathode 


surface. The erosion mechanism for the type I spot is beHeved 
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to be due to thermal evaporation. 


B. THEORY OF STATIONARY CATHODE SPOTS 

The theoretical modeling of cathode spots has progressed 
over the last 100 years, but up until the present time there 
has been no satisfactory understanding of the processes 
occurring within the spot. The difficulty occurs due to the 
diverse phenomena, which in themselves are extremely complicated 
to describe mathematically, that take place within the small 
cathode-plasma region known as a cathode spot. A _ detailed 
analysis would require an extensive use of the kinetic theories 
of Plasma Physics, Solid State Physics, Physics of High 
Temperature Phenomena, and Gas Dynamics. In order to avoid 
some of the complications the majority of models developed for 
the cathode spot are stationary, and therefore are unable to 
explain the initiation, motion, and disappearance of the cathode 
spot [Refs. 12 and 15]. 

The energy balance method developed by Lee and Greenwood 
fRef. 16] has been the basis of the majority of stationary 
models including those proposed by Kulyapin (Ref. 17], Kubono 
(Ref. 18], Beilis [{Ref. 19], and Moizhes and Nemchinskii [Refs. 20, 
21, 22, and 23). The latter work by Moizhes and Nemchinskii are 
the basis for the model developed in this paper. The cathode 


spot is divided into three regions as illustrated in Figure 4.1. 
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Figure 4.1 Model of the Cathode-Plasma Region 


Region I, referred to as the Langmuir region, is a thin 
Space charge sheath over which the cathode potential drop, Vk, 
occurs. It is in this region that the electrons emitted by the 
cathode are accelerated into the plasma and where the ions 
produced in the generation region are accelerated towards the 
cathode. Region II, referred to as the generation region or 
Knudsen Region, is where the neutral atoms evaporated from the 
cathode are ionized. Region HI is an expanding channel or 
transition region over which the highly dense current channel in 


regions I and II expands out into the bulk plasma. 
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L Fundamental Equations 

In developing a theoretical model of a cathode spot, one 
of the first decisions that must be made is the electron 
emission mechanism. In general most models use a combination 
of thermal and field equations that are expressed as 

Je=~f(T,-F) (eqn. 4.1 ) 
where je is the electron current density, Tk is the cathode 
temperature, and F is the electric field at the cathode surface. 
In actual use, expressions such as equation 2.41, 2.45, or 2.56 
are used. The particular equation used depends on the 
applicability under the given conditions. The electric field at 
the cathode surface is determined using MacKeown’s Equation, 
equation 2.36. 

An energy balance must exist in the plasma region above 


the cathode spot and can be described by equation 4.2. 


kT, 
e 


( eqn. 4.2 ) 





(eV ,+2kT,)= L(g ,+2kT,)+B 
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Energy is brought into the region by the electrons emitted from 
the cathode and consists of the thermal energy and the energy 
gained as the electrons were accelerated through the cathode 
potential drop. Energy is removed from the region by ions 
flowing back to the cathode. This energy loss consists of the 
ionization energy and the thermal energy of the ions. A further 
loss occurs from both ions and electrons that leave the region 


and pass into the bulk plasma and is described by the Peltier 
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heat flux /fkT,. In a plasma, in which Coulomb collisions are the 
predominant eneanine process, the energy transfer coefficient, 
8, is approximately equal to 3.21 It should be noted that 
electrons in the plasma, also, can leave the region towards the 
cathode, but only if energetic enough to over come the cathode 
potential drop, Vk. This reverse electron current can be 


approximated by 


| enaP | (eqn. 4.3 ) 
ae — | 
V2nkT , Kili 


For low plasma temperatures and high cathode potential drops 
the reverse electron current density is much less than the ion 
current density and can therefore be neglected in this modeL 

An energy balance at the cathode surface must also exist 


and can be described by equation 4.4. 
AK RT, = REE 2, + 0V,~6)- H(6+ 247 .)- gE ap | 
( eqn. 4.4 ) 
The predominant loss mechanisms considered are heat conduction, 
emissive cooling, and evaporation. The primary energy source 
considered is ion bombardment of the cathode surface. The ions 
bring to the surface thermal energy, ionization energy, and 
energy gained while being accelerated by the cathode potential 
drop. Energy is released, approximately equal to the work 
function, when ions and electrons recombine at the surface. The 


electrons emitted by the surface remove thermal energy and an 


45 


energy approximately equal to the work function. In addition, 
neutral ares evaporating from the surface, remove energy 
equivalent to the heat of vaporization. Heat conduction from 
the cathode spot is based on a circular spot of radius R at a 
constant temperature Tk and where K is the thermal conductivity 
of the cathode. The work function to be used in equation 4.4 is 


corrected for the Schottky effect and is given by 


a 7F 





= ¢,- eqn. 4.9 
=, ee ( eq ) 
The total current flow in the arc is given by 
f=nR*(j,+J;). (eqn. 4.6 ) 


A fraction of the neutral atoms evaporated from the cathode 
will return to the cathode. The return coefficient can be 
determined from equation 4.7 which is based on a Monte Carlo 


calculation performed by Nemchinskii [Ref. 24]. 


[SxT [Sa\" 
n= (Bea aaTt)( Beas =f : ( eqn. 4.6 ) 
p 


The ion current density is determined from 





I ae niit ( eqn. 4.8 ) 
where 
Pan 
- eqn. 4.9 
Py loa ( q ) 
and 
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5 /8kT. 
rade @eqn. 4.10 ) 


Substituting equations 4.9 and 4.10 into equation 4.8 gives the 


1on current density as 


aE 
eT eae 4a 
PR (Ta ( eq ) 


where P is the vapor pressure of the cathode material and «@ is 
the degree of ionization. The degree of ionization is determined 
from Saha’s Equation which is expressed as 


| 


a= ] + = aa “ ( =| Asal 
Peher, | eee Tey EID oithh 2 ) 


The erosion rate of neutral atoms is given by 
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ee (eqn. 4.13 
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The potential drop inside the quaSineutral plasma of the 


transition region is given by 


U=| (seen) Z)c *2)-BKT, ||] aa +2] 


( eqn. 4.14 ) 





where z is the average charge of the ions in the region, j is 


the total current density given by 


j=hitde C eqn. 4.15 ) 


and ¢ is the electric conductivity given by 
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7 Ts ( eqn. 4.16 ) 


where the Coulomb Logarithm is approximated as 


In A= 29 In( 2 *}. ( eqn. 4.17 ) 


p 
The total potential drop across the arc is given by 

Yor pe. (eqn. 4.18 ) 

2. Method of Solution 

The system of equations presented in the previous section 
is not closed. Therefore, more information is required in order 
to solve the system. This will be accomplished here using the 
Steenbeck Minimum Principle (Refs. 25 and 26}. Simply stated this 
requires that the total potential drop across the arc must be a 
minimum for a stationary arc to exist. Given an arc current and 
cathode temperature the system of equations is solved for the 
cathode potential drop, plasma temperature, ion current density, 
electron current density, erosion rate, potential drop in the 
transition region, and total arc potential drop. The cathode 
temperature that defines the stationary cathode spot, for the 
given arc current, is that cathode temperature that minimizes 
the total potential drop. If a minimum total arc potential can 
not be found for a given arc current, then a stable arc can not 


exist for this current. 
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The first step in the solution is to determine the plasma 
temperature, Tp, and cathode drop, Vk. To accomplish this a 
guess is made for the plasma temperature. The partial pressure 


of the cathode material is determined using [Ref. 27] 


A 
logP = ——*+ Blog? ,+CT,+D ( eqn. 4.19 ) 


PD 


mere i fommecoppeneA- —L/52x104,09 B= =24, C=40) and D= 13.21 ~The 
degree of ionization is determined uSing equation 4.12 and the 
return coefficient is determined using equation 4.6. The ion 
current density can now be determined using equation 4.1L 
Experimentally, it has been observed that the cathode drop is 
aperbcsamacely 15 volts, so this will be used as a first 
approximation in order to determine the electric field at the 
cathode surface using equation 2.36. The Schottky corrected 
work function is determined using equation 4.5 and the electron 
current density is determined using equation 4.1 The spot 
radius can now be determined using equation 4.6. The erosion 
rate is determined from equation 4.13. Equation 4.2 and 4.4 are 


each solved separately for the cathode potential which gives 


ii ' 
VaR “2 Ae, 2kT ,)* BRT (1 -4)-onr,} 
( eqn. 4.20 ) 


and 
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4 4 | ue g 
VY ..=—- T —({ — 2kT ,J-£&,-2kT —f : 
Ka PR], at ile JF proped ‘on 


( eqn. 4.21 ) 
The potential is then determined uSing both equation 4.20 and 
4.21. The results are compared and the plasma temperature is 
then adjusted. This procedure is then repeated until the 
potentials determined from both equations are in agreement. The 
conductivity is determined using equation 4.16 and 4.17. The 
potential drop in the transition region is then determined using 
equation 4.14. Finally, the total arc potential is determined 
from equation 4.18. Appendix J is a Fortran program, which uses 
the method described above, to calculate the cathode spot 


parameters. 
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V. COMPARTSON OF THEORY AND EXPERIMENTAL DATA 


The Fortran program listed in Appendix J was used to 
calculate the parameters of the cathode spot for arc currents 
of 100, 150, 175, 200, 300 and 400 Amperes. The program 
required knowledge of several material parameters which are 


specified in Table 1. 


TABLE 1 
CATHODE SPOT MODEL CONSTANTS 








Thermal Conductivity (W/m-kK) 


Figure 5.1 illustrates the variation of the total arc 


potential with cathode temperature for different arc currents. 
The cathode temperature for a stable cathode spot, and thus 
the parameters of the cathode spot, is determined from the 
minimum of the total arc potential It can be seen from Figure 
5.1 that a stable arc does not occur for are currents below 200 


Amperes. 
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Figure 5.1 Cathode Spot Stability 
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The results of the calculations are listed in Table 2 through 
7 for arc currents of 100, 150, 175, 200, 300, and 400 Amperes 
respectively. In addition Figures 5.2 through 5.7 illustrate the 
variation of the cathode potential, transition region potential 
and total arc potential for each arc current. 

The cathode spot parameters for stable arc currents of 200, 
300, and 400 Amperes are summarized in Table 8. The results 
are in general, order of magnitude, agreement with experimental 
data [{Ref. 23], however stability has been indicated at arc 
currents between 100 and 200 Amperes. This discrepancy occurs 
Since the material constants listed in Table 1 are not really 
constants at these temperatures. A more detailed comparison is 
not warranted, since the published experimental data for the 
cathode spot parameters vary widely. Confirmation of the models 
of the cathode phenomena is greatly limited by this lack of 


consistent experimental results. 


O3 













TABLE 2 
CATHODE SPOT PARAMETERS FOR A 100 AMPERE ARC ON A 
COPPER CATHODE 


Cathode 
Potential 
(Volts) 














Total Arc 
Potential 
(Volts) 


Transition 
Potential 
(Volts) 














Partial 
Pressure 
(atm.) 


Degree of 
Ionization 






Coefficient Radius 
(1072 sme 


Temperature 
(Kelvin) 
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TABLE 2 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 100 AMPERE ARC ON A 
COPPER CATHODE 


Cathode Electric Electron Ion Total 
Temperature Field Current Current Ciaran t 
(Kelvin) moe V7m) Density Density Density 
(102 A/m=) (102 A/m=) 108 A/m=) 


Cathode Ion/Total Coulomb eondtic— 
Temperature Current Logarithm Ca ey 
(Kelvin) Density (mho/m) 
Ratio 





wD 


POTENTIALS FOR A 100 AMPERE ARC 
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Figure 5.2 Potentials for a 100 Ampere Arc 
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TABLE 3 
CATHODE SPOT PARAMETERS FOR A 150 AMPERE ARC ON A 
VOPPEROCALHODE 


Cathode Plasma Cathode Transition Total Arc 
Temperature | Temperature Potential Potential Potential 
(Kelvin) (eV) (Volts) (Volts) (Volts) 


Cathode Partial Degree of Return Sipe 
Temperature Pressure Ionization Coefficient Radius 
(Kelvin) (atm.) CG 2 31) 


0.4716 
0.4742 
0.4768 
0.4794 
0.48721 
0.4849 
0.4877 
0.4907 
0.4938 
0.4971 
G.s0Cs 
0.5044 
0.5084 
Ceo 7 
Grane 
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TABLE 3 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 150 AMPERE ARC ON A 
COPPER CATHODE 
















Ion Total 


Cathode Electric Electron 







Temperature Field Current Current Current 
(Kelvin) (lO V/m) Density Density Density 
(LOS A/m=) (102 A/m=) (108 A/m=) 





Ca ceeds [en eead Coulomb Senne — Erosion 
Temperature Current Logarithm tivity Rate 
(Kelvin) Density (mho/m) TOo=24 
Ratio atoms/m<-—s) 


0.5443 
0.5194 
0.4946 
0.4700 
0.4458 
0.4221 
ORS Shy, 
0.3762 
0.3542 
OLSS29 
O.o1Z2 
OeZz7 2S 
On 2ele2 2 
0.2542 
OrzSor 
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FOTENTIALS FOR A 150 
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Potentials for a 150 Ampere Arc 


Figure 5.3 
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TABLE 4 


CATHODE SPOT PARAMETERS FOR A 175 AMPERE ARC ON A 
COPPER CATHODE 


Cathode Transition Total Arc 
Potential Potential Potential 
(Volts) (Volts) (Volts) 





























Cathode Plasma 
Temperature |]Temperature 
(Kelvin) (eV) 








eo 
Ze 57, 
135.04 
nS. 7 2 
14.39 
15.04 
15.66 
16.24 
oe? o 
M21 
oS 
7.6 
17 .o0 
V7 6S 
l/ais 















Return 
Coefficient 


Cathode Partial 


Degree of 
Ionization 


Spot 
Radius 
(Om? m 


Temperature Pressure 





(Kelvin) (atm.) 





O-ls7e 
Onse2 
Oz750 
0.1842 
Ono Ss 
0.2040 
0.2148 
0.2264 


O22 597 
Oe a25 
Or2o76 
0.2843 
O2SOS! 
0.3246 
0.3495 
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0.4714 
0.4739 
0.4765 
0.4791 
0.4818 
0.4845 
0.4874 
0.4903 
0.4934 
0.4967 
O25002 
Ons 57 
0.3080 
Orsi 5 
0.3175 




















TABLE 4 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 175 AMPERE ARC ON A 
COPPER CATHODE 

















Cathode 










Electric Electron Ion otal 



























Temperature Field Current Current Current 
(Kelvin) (LOS V/m) Density Density Density 
(108 A/m=) (108% A/m=) (102 A/m=) 


Schade Tenarote | Boutiont Gerduet Erosion 
Temperature Current Logarithm tivity Rate 
(Kelvin) Density (mho/m) Om 
Ratio atoms/m=-s) 


0.5448 
O-a sc 
0.4940 
0.4695 
0.4453 
0.4216 
0.3984 
OnS 755 
OS 57 
O.3S26 
CrS1Z0 
OZ 721 
O27 Ze 
0.2542 
O22 561 


orl 


POTENTIALS FOR A 175 AMPERE ARC 
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Potentials for a 175 Ampere Arc 


Figure 5.4 
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TABLE 5 
CATHODE SPOT PARAMETERS FOR A 200 AMPERE ARC ON A 
COPPER CATHODE 


Cathode Plasma Cathode Transition Total Arc 
Temperature | Temperature Potential Potential Potential 
(Kelvin) (eV) (Volts) (Volts) (Volts) 


Cathode Partial Degree of Return Spot 
Temperature Pressure Tonization Coefficient Radius 
(Kelvin) (atm.) Om m,) 


0.4713 
0.4736 
0.4762 
0.4788 
0.4815 
0.4842 
0.4871 
0.4900 
Oea7 Si 
0.4964 
0.4998 
O-50S6 
OfD076 
Ora iz | 
One T7il 
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TABLE 5 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 200 AMPERE ARC ON A 
COPPER CATHODE 


Electron 
























Total 







Cathode Electric Ion 









Temperature Field Current Current Current 
(Kelvin) (LOS V/m) Density Density Density 
(Oe A/m=) 102 A/m=) OSA m=) 


Eamode Ion/Total Coulomb Conduc— Erosion 
Temperature Current Logarithm aw ney Rate 


(Kelvin) Density (mho/m) (1O=7 
Ratio atoms/m=-s) 


0.9438 
Ozelss 
0.4935 
0.4690 
0.4449 
0.4212 
0.3981 
O25 755 
0.3536 
Ono 52S 
0.3118 
Orzo 17 
O.27 2 
0.2541 
Osz2 560 
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POTENTIALS FOR A 200 AMPERE ARC 
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Figure 5.5 Potentials for a 200 Ampere Arc 
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: TABLE 6 
CATHODE SPOT PARAMETERS FOR A 300 AMPERE ARC ON A 
COPPER CATHODE 


Cathode 
Potential 
(Volts) 


Total Arc 
Potential 
(Volts) 


Transition 
Potential 
(Volts) 


Cathode 
Temperature | Temperature 
(Kelvin) 








Zee 
21.06 
19270 
18.86 
17.94 
17 hs 
16.41 
Laer oS 
15.24 
14.77 
14.38 
14.07 
13.82 
US {ops 
Nees >) 




















Degree of 
Tonization 


Cathode 
Temperature 
(Kelvin) 


Partial 
Pressure 
(atm.) 


Spot 





Radius 
(1OF>* “my 


Coefficient 








Ola 0 
©O2162Z6 
Oey 12 
0.1801 
GAS576 
Oac7S 
0.2101 
0.2214 


0.25336 
0.2470 
O2Z016 
Om 7 7 
OEZz ToS 
Ons 2 
0.3415 
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0.4706 
0.4730 
0.4755 
0.4781 
0.4808 
0.4835 
0.4863 
0.4892 
0.4923 
0.4955 
0.4990 
O20 27 
O25067 
eo i 
Oe LoO 








TABLE 6 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 300 AMPERE ARC ON A 
COPPER CATHODE 


Cathode Electric Electron Ion Total 
Temperature Field Current Current Current 
(Kelvin) (10S V/m) Density Density Density 
(1108S A/m=) (102 A/m=) 108 A/m=) 


Pode eet Eotiens Boe Eee 
Temperature Current Logarithm tAVity Rate 
(Kelvin) Density (mho/m) dO=¢@ 
Ratio atoms/m=—s) 


0.5423 
©7516 7 
0.4922 
0.4678 
0.4438 
0.4202 
Caso 71 
0.3747 
O27 > 
CES Si7 
C113 
C2271 
0.2724 
0.29357 
G-2 500 
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200 AMPERE ARC 


FOTENTIALS FOR A 
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Potentials for a 300 Ampere Arc 


Figure 5.6 
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TABLE 7 
RS FOR A 400 AMPERE ARC ON A 
COPPER CATHODE 


Cathode 
Potential 
(Volts) 

















Degree of 
Ionization 





Onlaso 
O.1610 
Oue7 > 
0.1784 
Onley 7 
Onl 7c 
O-2.0 a 
Onl S 


0.2314 
0.2446 
0.2571 
O.27 32 
0.2934 
0.3141 
0.3380 
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Transition 
Potential 


(Volts) 





Coefficient 


0.4702 
0.4726 
0.4751 
0.4777 
0.4803 
0.4830 
0.4858 
0.4888 
0.4918 
0.4950 
0.4984 
0.5021 
O.2001 
O.a105 
0.5154 


Total Arc 








TABLE 7 (CONTINUED) 
CATHODE SPOT PARAMETERS FOR A 400 AMPERE ARC ON A 
COPPER CATHODE 


Electron Ion Total 
Current Current Current 


Cathode Electric 
Temperature Field 
(Kelvin) (1lO8 V/m) 


SOOO 
SS a0 
3600 
3690 
3700 
S750 
3800 


3850 
3900 
S750 
4000 
4050 
4100 
4150 








2.740 
ceo 
3.456 
3.7435 
4.049 
4.377 
4.728 


S107 
Seale 
S701 
6.446 
oe7 oO 
Tea7 2 
8.234 





Density Density Density 
(108 A/m=) (102 A/m=) (108 A/m=) 
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TABLE 8 
CATHODE SPOT PARAMETERS FOR STABLE ARCS 


200 





Parameter 


Sen EELS 


3700 3600 
a eC 
lSne © 17.94 19.81 


ie 
3.022 2) ASG, Tez 
gaol? 1.514 1.036 
ao moo 


ee ee 


| 


ta 


Cathode Temperature 3850 


(Kelvin) 
Plasma Temperature (eV) 
Cathode Potential (Volts) 






Spot Radius (10-4 Mm) 


Ion Current Density (108 
A/m2 ) 


Rlectron Current Density 
(108 A/m2 ) 


Total Current Density (108 
A/m2 ) 


Brosion Rate (1027 
atoms /m2 -s) 
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fr -4.072 
6.970 412 2.108 
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VI. CONCLUSTONS AND RECOMMENDATIONS 


This study has attempted to answer some of the questions 
associated with the cathode phenomena in a unipolar arc. 
However, a more detailed analysis of the cathode spot must be 
accomplished before a consistent understanding of the diverse 
phenomena occurring within the arc can be achieved. The 
stationary model presented her is only a beginning and a tool 
for the development of future models. It can not answer the 
questions regarding arc initiation, extinction, and spot motion 
that will have to be addressed in future models if they are to 
be accepted. Additional research with stationary models should 
consider including the reverse electron current from the plasma 
back to the metal changing from a constant cathode temperature 
in the arc radius to a constant heat flux with a varying 
temperature, considering the ohmic heating of the cathode spot 
by the electron current, and varying the cathode material 
parameters to determine their influence on the arc stability. It 
is hoped that this study will help to stimulate further interest 


and research into the cathode spot of a unipolar arc. 
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APPENDIX A GENERAL EMISSION PROGRAM 


FOCI KOKORO OOK OOO OOK OKO KOK KKK KOK KOK KKK KOK 


« EMM1. EXE 
x 

‘ PURPOSE 

. CALCULATE THE ELECTRON CURRENT DENSITY 

so EMITTED FOR A RANGE OF ELECTRIC FIELDS AT 
c THE SURFACE OF A CATHODE WITH A GIVEN WORK 
so FUNCTION AND CATHODE TEMPERATURE. 

*K 

3 DESCRIPTION OF VARIABLES 

i; CATHODE TEMPERATURE 

F ELECTRIC FIELD 

FO INITIAL ELECTRIC FIELD 

. PHI WORK FUNCTION 

« DF ELECTRIC FIELD INCREMENT 

* J CURRENT DENSITY 

. ECODE ERROR CODE 

. 0 = GOOD 

. 1 = ERROR 

ss TYP TYPE OF EMISSION 

* i= THERMIONIC 

. 2= INTERMEDIATE 

. 3= FIELD 

. = NONE 

K 

SUBROUTINES REQUIRED 

. EMMAPP(T, F, PHI, TYP) 

TFEMMUT, F, PHI, J, ECODE) 

. IF TEMM1(T, F, PHI, J, ECODE) 

. FTEMMUT,F, PHI, J, ECODE) 

* 

KOKOK OK KK KK KK KKK KOK KOK KOK KKK OK KK KKK KK OK KK KK KK OK KOK OK OK OK KK OK OK KOK 


REAL*8 T,F,FO,PHI,DF,J 

INTEGER I,N,ECODE,TYP 

OPEN(4, FILE= CON’) 
OPEN(6,FILE="EMM1LDAT’,STATUS="NEW’) 

WRITE(*,*) “ENTER THE NOMBER OF DATA POINTS: © 
READ(*,*) N 

WRITE(*,*) “ENTER THE ABSOLUTE TEMPERATURE (K): ° 
READ(*,*) T 

WRITE(*,*) “ENTER INITIAL ELECTRIC FIELD (V/M): ° 
READ(*,*) FO 

WRITE(*,*) ENTER ELECTRIC FIELD INCREMENT (V/M): ~ 
READ(*,*) DF 


aa 


10 
20) 


900 


WRITE(*,*) “ENTER WORK FUNCTION (eV): ° 

READ(*,*) PHI 

F=FO 

J=0. 

DO 20 I=1,N 
CALL EMMAPP(T, FS PHI? 72) 
IF(TYP.EQ.1) CALL TREMMIG TE Foi, J, BeOpR) 
IF(TYP.EQ.2) CALL IF TEMMIGE SE PHig 2eoun) 
TE(TYP.EQ.3) CALL ETEMMIGT TE PHI, 2 CODE) 
IF(TYP.EQ@.4) GOTO 10 
WRITE(6,900) F, J, ECORE wie 
F=F + DF 

CONTINUE 

a1 OF 

FORMAT(1X, D12.6,5X,D12.6,5X,15,5X, I5) 

END 
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APPENDIX B GENERAL EMISSION SUBROUTINE 
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* 
*K 
x 
* 
* 
* 
K 
* 
*K 
* 
* 
x 
* 
* 
* 
* 
x 
* 
x 
2K 
K 
Ps 
* 
i 
*K 
* 
K 
* 
x 


SUBROUTINE EMM 


PURPOSE 


CALCULATE THE ELECTRON CURRENT DENSITY 

FOR A GIVEN ELECTRIC FIELD AT THE SURFACE 
OF A CATHODE WITH A GIVEN WORK FUNCTION AND 
CATHODE TEMPERATURE. 


DESCRIPTION OF VARIABLES 


ls 

F 

Pl 

J 
ECODE 


Pe 


CATHODE TEMPERATURE 
ELECTRIC FIELD 

WORK FUNCTION 
CURRENT DENSITY 


ERROR CODE 
0 = GOOD 
ee ee ROR 


TY PESO less LON 
l= THERMIONIC 
2= INTERMEDIATE 
3= FIELD 
4= NONE 


SUBROUTINES REQUIRED 
Hive leva Meals soUladg aay) 
MEM Bb Pal, CODE) 
IFTEMMU(T,F, PHI,J, ECODE) 
FTEMMiU(T,F, PHI,J,ECODE) 


KOKO OOOO OOK OOK IOI OK OKO KOK IOIOK FOI OK OOK OK KOK HOOK OK OOK 


SWEROULINE EMM(T,', PHI,J,ECODE, TYP) 
REAL*8 T,F,PHI,J 
PN EGGeR ECODE, TYP 

CAPE EMMAPP(T,F,PHI,TYP) 

MeCiyeehQ 1) CALL TFEMMI(T,F, PHI, RCODE) 
IF(TYP.EQ.2) CALL IFTEMMI(T,F, PHI,J, ECODE) 
IF(TYP.EQ.3) CALL FTEMMIU(T,F, PHI,J, ECODE) 
RETURN 

END 


co 
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APPENDIX C EMISSION APPLICABILITY SUBROUTINE 
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OROOIOK OK OKI IOOK OOOO OOO IOI OKO OOOO KOK KOK OK KOK KOK 


SUBROUTINES REQUIRED 
NONE 


x 

* SUBROUTINE EMMAPP 

* 

x PURPOSE 

* DETERMINES APPLICABILITY OF THE SUBROUTINES 
r TFEMM1, IFTEMM1, AND FTEMM1 

x 

x DESCRIPTION OF VARIABLES 

* T CATHODE TEMPERATURE 

* FO ELECTRIC FIELD AT CATHODE SURFACE 
x PHIO WORK FUNCTION 

x TYP TYPE OF EMISSION 

* 1= THERMIONIC 

* 2= INTERMEDIATE 

* 3= FIELD 

x 4= NONE 

x 

x 

x 

> 4 


ROR IOKOKOKOKOKOK OOK KOK OK OKOOROIOROK KOK OK KOK OK OK OK OK KOK OK OK OK OK KKK OK 


SUBROUTINE EMMAPP(T bOr ai eas at) 
REAL*8 T,f,PHI,M,E,K, Hj Bogs eo Use eee 
REAL*8 FO,E1,F1,KTO\ Kili) eee Een, BT AO 
REAL*8 TEST1,TEST2,TEST3,TEST4 
INTEGER ECODE1,ECODEZ2,ECODE3,TYP,N 
M= Judo Eaot 

E=16022E-19 

K=1.3807E-23 

H=1.0546E-34 

EPS=8.8542E-12 

PI=3.141592654 

E1=(4* PIKE PS*H)*k*2 /(MKEX* 4) 
F1=K1**2 *E**3 /(4*PTKEPS) 

P=RO*R1 

KTO=K*T 

KT=KTO*E1 

PHI=PHIO *E1*E 

N=0 

D=F**(3./4.)/(PI*XKT) 

PG Daur aun le 

TEST1=DLOG((1- D)/D)-1/(D*(1- D)) 
TEST2=- PI/(E**(3. /4.))*(PHI- DSQRT(F)) 
THES 13s DEOG(I— ») 7) 7 Ge) 
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20 
50 


40 


50 


TEST4=-PI/(F*x(1. /8.)) 
IF((TEST1GT.TEST2). AND.(TEST3.GT.TEST4)) THEN 
ECODE1=0 
ELSE 
Peo DR I= 1 
ENDIF 
SOO 30 
ECODE1=1 
ETAOQ=-F**2 /((8*KT**2) 
Y=DSQRT(F)/PHI 
CAli YY Gi val Yoke 1) 
YE=DSQRT(F)/(-ETAO) 
OA DV Gan Bil BER ) 
ETA=-FR*2 /(8*KT**2 *TE**2) 
IF(DABS(ETA-ETAO).GT. 10D-6) THEN 
KTAOSETA-(ETA-ETAO)/2.1 
N=N+1 
LP (NIG 200)) THEN 
HCO Z=1 
GOsee 50 


D=2*DSQRT(2/YE)*TE/PI 
TESTI=1/YE 
esa tebe x x1 /4.)*D/(P1*(D-1)) 
TEST3=-F**2 /(8XKT**K2*TE*K*Z) 
TEST 4=- PHI+KT A1- F/A2* DSQRT(2* PHI)*KT *T Y)) 
ito alt hsl2) AND. (TEST3.GI-TEST4)) THEN 
ECODE2=0 
ELSE 
ECODE2=1 
ENDIF 
D=(2* DSQRT (2.0 * PHI)*PI*KT)/F 
Di=4*DSQRT(2 *PHI)*PHI/(3*F) 
TEST1=PHI- DSQRT(F) 
TEST2Z=Fx*x*(3./4.)/PI + KT /(i- D*TY/PD 
iio =e YP 1 
TEST4=KT*DSQRT((V* PHI/F)* DSQRT(2 * PHD) /(PH1**2 - F)) 
PEO EEG! TEStTZ)AND.TEST3.GT.TEST4)) THEN 
ECODE3=0 
ELSE 
ECODE3=1 
ENDIF 
IF(ECODELEQ.0) THEN 
iy Ea 
ELSE 
IF(ECODE2.EQ.0) THEN 
tera 2 
ELSE 
IF(ECODE3.EQ.0) THEN 
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RETURN 
END 
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APPENDIX D T-F BMISSION SUBROUTINE 
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x 

x SUBROUTINE TFEMM1 

*K 

* PURE Cos 

* CALCULATES THE ELECTRON CURRENT DENSITY 
* FOR THERMIONIC FIELD EMISSION 

4 

* DESCRIPTION OF VARIABLES 

* at CATHODE TEMPERATURE 

* FO ELECTRIC FIELD AT CATHODE SURFACE 
* Eneae WORK FUNCTION 

*K JO CURRENT DENSITY 

x ECODE ERROR Gorn 

* 0= GOor 

*K l= ERROR 

. 

x SUBROUTINES REQUIRED 

x NONE 

* 


OIOIOIOIOIOIOIOIOKIOICIOIOICIOICIOFRCICIOIO I IOC IK OOOO OOK OOOO OOK OK ACK CK 


SUBROUTINE TFEMMi(T,FO, PHIO,J0, ECODE) 
REAL*8 T,F,PHI,J,M,E,K,H,EPS, PI, PHIO,D 
REAL*8 FO,J0,E1,F1i,J1,KTO,KT 
Reao*«s THoll,TEST2,TEST3,TEST4 
INTEGER ECODE 
M=9.1095E-31 
E=1.6022E-19 
K=1.3807E-23 
H=1.0546E-34 
EPS=8.8542E-12 
El=3-141592654 
Bi=(4*PIXEPS *H)*k*2 /(MKE**4) 
Fi=E1**2*E**3 /(4*PIXEPS) 
J1i=E1**2 *H**3 /(M*E) 
F=FOXF1 
Kl 02 KT 
KT=KTO*EI1 
PHI=PHIO *E1xE 
D=F*x«(3./4.)/(PI*KT) 
J=(KT**2 (2 * PI* «2 ))x((PI* D)/DS IN(PI*D))* 
1 DEXP(-(PHI-DSQRT(F))/KT) 
JIGS Ii 
IF(i-D) 20,20,10 
10 TEST1i=DLOG((1- D)/D)-1/(D*(1- D)) 
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20 


TEST2=-PIK(F* (3. /4.))*(PHI- DSQRT(F)) 
TEST3=DLOG((1i- D)/D)- 1/(1- D) 
TEST4=-PI/(E*x(1. /8.)) 
IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4)) THEN 
ECODE=0 
ELSE 
ECODE=1 
ENDIF 
RETURN 
ECODE=1 
RETURN 
END 
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APPENDIX KE FIKLD-THERMIONIC EMISSION SUBROUTINE 


ROCIO OOOO KOKORO OOK OOK KOK OK KOK KK KOK KOK KOK OK OK OK OK OK OK OK KOK KOK KOK 


* SUBROUTINE FTEMM1 

K 

* PURPOSE 

* CALCULATES THE ELECTRON CURRENT DENSITY 

x FOR FIELD-THERMIONIC EMISSION 

«x 

re DESCRIPTION OF VARIABLES 

* iT CATHODE TEMPERATURE 

* FO ELECTRIC FIELD AT CATHODE SURFACE 
* PHIO WORK FUNCTION 

« JO CURRENT DENSITY 

* ECODE ERROR CODE 

* O= GOOD 

* l= ERROR 

«K 

* SUBROUTINES REQUIRED 

é VTY(Y, V, TY, IER) 

* 

OK OK KOK OK KOKO OK OK OK OK KOK KOK KOK OK KKK OK OK OK OK OK OK OK KOK OK OK OK OK OK OK KOK OK OK OK OK OK KK OK OK OK 


SUBROUTINE FTEMMI(T,FO,PHI0O,J0,ECODE) 
Po Omens, tah, H bro, bl, PRI0eD bi, TY 
Reb omeoll, | KolZ, TEST 3; PeSb4ayy, ¥,61,62,C3 
Ree beer Oj 0 ,H 1, Pilwl,KTO,KT 

INTEGER ECODE 

ey Jono 1. 

B= s0228-19 

K=L3807E-23 

H=10546E-34 

Hro=8.S5420-12 

PI=3.141592654 

H1=(4* PIKEPS*H)**2 /(M*KE**4) 
PI=HI*X*K2*E**K3 /(4*PIXEPS) 

J1I=H1**2Z XH**3 (ME) 

F=FQOX*F1 

KTO=K*T 

KT=KTO*E1 

PHI=PHIO *E1*E 
D=(2*DSQRT(2.0 * PHI)* PIXKT)/F 
D1l=4*DSQRT(2* PHI)* PHI/(3 *F) 
Y=DSQRT(F)/PHI 

CALL VIY(Y,V,TY,IER) 

J=E**2 *(D*XTY/DSIN(D*TY))* DEXP(- D1*V)/(16 * PIk*2 * PHIXT Y**2 ) 
JO=JS/J1 
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TEST1=PHI- DSQRT(F) 
TEST2=F**(3./4.)/PI + KT/(1-D*TY/PI) 
TEST3=1-D*TY/PI 
TEST4=KT*DSQRT((V*PHI/F)* DSQRT(2* PHI) /(PHI*2 - F)) 
IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4)) THEN 
ECODE=0 
ELSE 
ECODE=1 
ENDIF 
RETURN 
END 
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APPENDIX F COMPLETE ELLIPTIC INTEGRAL OF THE 1ST 
KIND 


OOOO COCO CK OOOO IK OK KK 
SUBROUTINE CEL1 
PURPOSE 
CALCULATE COMPLETE ELLIPTIC INTEGRAL OF FIRST 
KIND 


USAGE 
CALL CEL1(RES, AK, IER) 


DESCRIPTION OF PARAMETERS 


RES = RHOULT VALUE 
AK = HOU EbUS “INEGH 
IER - RESULTANT ERROR CODE WHERE 


IER=0 NO ERROR 
IER=1 AK NOT IN RANGE -1 TO +1 


REMARKS 
BoB RESULT IS SET TO LE7S IF ABS(AK) GE 1 
FOR MODULUS AK AND COMPLEMENTARY MODULUS CK, 
EQUATION AK*AK+CK*CK=10 IS USED. 
AK MUST BE IN THE RANGE -1 TO #+1 


SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 


METHOD 
DEFINITION 
CEL1(AK)= INTEGRAL(1/SQRT((1+T *T)*(1+(CKxXT)**2)), 
SUMMED OVER T FROM 0 TO INFINITY). 
EQUIVALENT ARE THE DEFINITIONS 
CELI(AK)=INTEGRAL(1/((COS(T)SQRT(1+(CK*TAN(T))**2)) 
SME peOvER T FROM © TO Pi72), 
CEL1(AK)=INTEGRAL(1/SQRT(1-(AK*SIN(T))**2), 
Sonne OVERS) ROM Oe TO Piz), WHERE 
K=SQRT(1.-CK*CK). 
EVALUATION 
LANDENS TRANSFORMATION IS USED FOR CALCULATION. 
REFERENCE 
R.BULIRSCH, “NUMERICAL CALCULATION OF ELLIPTIC 
INTEGRALS AND ELLIPTIC FUNCTIONS’, HANDBOOK SERIES 
SPECIAL FUNCTIONS, NUMERISCHE MATHEMATIK VOL. 7, 
Ios ch. 1O- 90. 


Cy 


= C22 


KOI OOK IOI IOI OK OKO OK IK IOI ICICI IOI OK ICO’ OOK OKO OK IK KK 2K 


SUBROUTINE CELI1(RES, AK, IER) 
REAL*8 ARI,GEO,AK, RES, AARI 
IER=0 

ARI=2. 

GEO=(0.5- AK)+0.5 
GEO=GEO+GEOxAK 
RES=0.5 

IF(GEO)1,2, 4 

TER=1 

RES=LE38 

RETURN 

GEO=GEOX*AARI 
GEO=DSQRT(GEO) 
GEO=GEO0+GEO 

AARI=ARI 

ARI=ARI+GEO 
RES=RES+RES 
IF(GEO/ABRI-0.9999)35S%5 
RES=RES /ARI*6.283185E0 
RETURN 

END 
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APPENDIX G COMPLETE ELLIPTIC INTEGRAL OF 2ND KIND 


ORCC OOK OOOO OOK OOK OOK OK 2K KK KK 
SUBROUTINE CEL2 
PURPOSE 
COMPUTES THE GENERALIZED COMPLETE ELLIPTIC 
INTEGRAL OF SECOND KIND. 


USAGE 
CALL CEL2(RES, AK, A,B, IER) 


DESCRIPTION OF PARAMETERS 


RES - RESULT VALUE 

AK ——sOpULUS (INEUR) 

A - CONSTANT TERM IN NUMERATOR 

B - FACTOR OF QUADRATIC TERM IN NUMERATOR 
IER - RESULTANT ERROR CODE WHERE 


TER=0 NO ERROR 
1HR=1° AK NOT IN RANGES T.O +] 


REMARKS 
POR AbotAt GR 1 THE RESULT IS SHIShO LE/S IF B 
Ie Oomev ETO =LEf> Ih BIS NEGATIVE: 
SPECIAL CASES ARE 
K(K) OBTAINED WITH A = 
E(K) OBTAINED WITH A = 
COMPLEMENTARY MODULUS. 
B(K) OBTAINED WITH A = 1, 0 
D(K) OBTAINED WITH A = 0, B 1 
Wheeceene fy 6, D DEN INES SPECIAL CASES Ob THE 
GENERALIZED COMPLETE ELLIPTIC INTEGRAL OF SECOND 
KIND IN THE USUAL NOTATION, AND THE ARGUMENT K OF 
THESE FUNCTIONS MEANS THE MODOLOUS. 


1, 
I, 


1 
CK*XCK WHERE CK IS 


wo Ww 


SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 


METHOD 
DEFINITION 
RES=INTEGRAL((A+B*T*T)ASQRT((1+T*T) 
Re Gier  AZ))SCe)) 

SUMMED OVER | ROM 0 TOsINE INIT Y): 

EVALUATION 

LANDENS TRANSFORMATION IS USED FOR CALCULATION. 
REFERENCE 
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R.BULIRSCH, “NUMERICAL CALCULATION OF ELLIPTIC 
INTEGRALS AND ELLIPTIC FUNCTIONS’, HANDBOOK SERIES 
SPECIAL FUNCTIONS, NUMERISCHE MATHEMATIK VOL. 7, 
ibeksis)) deles fteles)0. 


ROOK OO KOO IOIOICOIOIOICICICIOIOR OO OOOO OK OKO OOK OK OOK KOK 


SUBROUTINE CEL2(RES, AK, A,B, IER) 
REAL*8 RES,AK,A,B, ARI, GEO, Al, BO, AARI 
IER=0 

ARI=2. 

GEO=(0.5-AK)+0.5 
GEO=GEO+GEOxAK 

RES=A 

A1=A+B 

BO=B+B 

IF(GEO)1,2,6 

IER=1 

IF(B)3,8,4 

RES=-1LE38 

RETURN 

RES=LE38 

RETURN 

GEO=GEOx* AARI 

GEO= DSQRT(GEO) 
GEO=GEO0+GEO 

AARI=ARI 

ARI= ARI+GEO 

BO =BO+RES*GEO 

RES=Al1 

BO = BO + BO 

Ai=BO /ARI+Al1 

IF (GEO7AARI-0) 99 92)5) ia 
RES=A1/ARI 
RES=RES+0.5707963E0*RES 
RETURN 

END 
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APPENDIX H FOWLER-NORDHIEM ELLIPTIC FUNCTIONS 


OCC OOK COCO IOI OK CK KOK 
SUBROUTINE VTY 


jel Ne let O}S)} 0 
CALCULATES THE FOWLER-NORDHIEM ELLIPTIC 
FUNCTIONS 


* 
*K 

x« 

*K 

K 

*K 

* DESCRIPTION OF VARIABLES 

* Y FUNCTION ARGUMENT 

V FOWLER-NORDHIEM ELLIPTIC V 
* TY FOWLER-NORDHIEM ELLIPTIC T 
* IER ERROR CODE 

* O= GOOD 

* 1= ERROR 
aK 

*K 

* 

*K 

*K 

«K 


SUBROUTINES REQUIRED 
CELIU(RES, AK, IER) 
CEL2(RES, AK, A, B, IER) 


ROOK OOK IO IO OOK OOK OOOO OR OK IK OK KOK OK OK OK OK OK OK KOK OK OK 


SUSROUTINE Vie Vy, TY, TER) 
Rime Soey Vo, C1, C2,C3 
HINGE Gi Repel Gye del BZ 
TE(Y.GT] 10) THEN 
CISRSQR hay —1)/(2*Y )) 
CALL CEL1(C2,C1,{ER1) 
CALL CEL2(C3,C1,10D00,10D00-C1*C1, IER2) 
V=-DSQRT(Y /2)*(-2*C3 + (Y+1)*CZ) 
ert /DSORTOUEY elt y)*xCa — Ye) 
ELSE 
C1=DSQRT((1- Y)A(1+ Y)) 
CAG CHEWCZ,C1, TERI) 
CAE CEEZ(C3,C1, 1 0D00, PGpigeeircl TeRZ) 
V=DSQRT(1+ Y)*(C3-Y*C2) 
imi Doakitir jet yjxC3 = Y*CZ) 
ENDIF 
IF ((TER1.EQ.1).0R.(IER2.EQ.1)) THEN 
IE R=1 


g1 
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APPENDIX I INTERMEDIATE F-T EMISSION SUBROUTINE 
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* SUBROUTINE IFTEMM1 
*K 

* PURPOSE 

* CALCULATES THE ELECTRON CURRENT DENSITY 

* FOR INTERMEDIATE FIELD-THERMIONIC EMISSION 
* 

* DESCRIPTION OF VARIABLES 

x T CATHODE TEMPERATURE 

* FO ELECTRIC FIELD AT CATHODE SURFACE 
x PHIO WORK FUNCTION 

* JO CURRENT DENSITY 

* ECODE ERROR CODE 

* O0= GOOD 

* 1= ERROR 

*K 

* SUBROUTINES REQUIRED 

* VTY(Y, V, TY, IER) 

x 

OOK OOK ICICI KOO OOOO OK OK OK KKK OK OK 


SUBROUTINE IFTEMMI(T,FO,PHIO,J0,ECODE) 
REAL*S Dt, PHI,J,M,E,K,H, EPS; Pie oreo 1 
REAL*8 TEST1,TEST2,TEST3, TEST4, Y,V,C1,C2,C3,VE 
REAL*8 FO,J0,E1,F1,J1,KTO,KT,ETA,ETAO,YE,THETA 
INTEGER ECODE 

Mao 109onS 3d 

BE=160 222-19 

Ke PSsO7h=23 

H=10546E-34 

EPS=8.8542E-12 

PI=3.141592654 

Bi=(4* PIKEPS*H)**2 /(M*KE**4) 
P1=H1x*2 *E**3 /(4*PIXEPS) 

J1=E1**2*H**3 /(M*KE) 

F=FOX*F1 

KTO=KxT 

KT=KTO*E1 

PHI=PHIO *E1*E 

ETAQ=-F**2 /(8 *KT**2) 

Y=DSQRT(F)/PHI 

CALI Vira ys ly tek 

YE=DSQRT(F)/(- ETAO ) 

CALL VTY(CYE, VE, TE, IER2) 

ETA=-EFR*2 /((8*KT**X2*TEXKXK2) 


SA 


IF(DABS(ETA-ETAO).GT. 10D-6) THEN 
ETAQ=ETA-(ETA-ETAO)/2.1 
GOTO 20 

RNDIF 

THETA=3 (TE**2) - 2*VE/(TE**3) 

J=F*DSQRT(KTATE A2*PI))/(2*PD*DEXP(-PHI/KT + F*x*2* 

1 THETA/(24*KT**3)) 

JO=JS/J1 

D=2*DSQRT(2/YE)*TE/PI 

eS ely i 

TEST2=1 + F**(1./4.)*D/(PI*(D-1)) 

TEST3=-FR*2 /((BXKTKXZ*KTEXK*2) 

TEST4=-— PHI+KT A1-F /(2* DSQRT(2* PHID*XKT*TY)) 

iit hol eels 2 yANDTEST3 Gil. beSt4)) THEN 
BCODE=0 

BLSE 
KCODE=1 

ENDIF 

RETURN 

END 
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APPENDIX J ARC PROGRAM 
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ARC. EXE 


PURPOSE 


DETERMINES THE CHARACTERISTIC VALUES FOR 
THE CATHODE SPOT OF AN ARC BURNING ON A 
COPPER PLATE 


DESCRIPTION OF VARIABLES 


ie 


er 
Eells 


SIG 
avis 


VK 
VT 


JE 
JI 
KTK 
Kal 


KKH HHKHHRHHKEHKEHHEHKHKEHHKKEKHKKKHKEHHE KKH HE HHKHKEHNHEHHHKHHRHH RHR HX 
NN 


HEAT OF VAPORIZATION 

PARTIAL PRESSURE CONSTANT 
MASS OF COPPER ATOM 

IONIZATION POTENTIAL ENERGY 
WORK FUNCTION (EV) 

FRACTION ENERGY REMOVED WITH 
COULOMB INTERACTION 

THERMAL CONDUCTIVITY 

ERROSION ENERGY 

AVERAGE ION CHARGE 

AVOGADRO’S NUMBER 

ELECTRON MASS 

ELECTRON CHARGE 

BOLTZMAN’S CONSTANT 

PLANCK’S CONSTANT 

PERMITIVITY OF FREE SPACE 
DEGREE OF IONIZATION 
VAPORIZATION RETURN COEFFICIENT 
ELECTRIC FIELD AT CATHODE SURFACE 
EROSION RATE 

SCHOTTKY CORRECTED WORK FUNCTION 
EE 

VAPOR PARTIAL PRESSURE 

RADIUS OF CATHODE SPOT 

PLASMA CONDUCTIVITY 

CATHODE TEMPERATURE (K) 
POTENTIAL DROP IN QUASINEUTRAL 
PLASMA 

CATHODE POTENTIAL DROP 

TOTAL POTENTIAL DROP ACROSS ARC 
TOTAL CURRENT DENSITY 
ELECTRON CURRENT DENSITY 

ION CURRENT DENSITY 

CATHODE TEMPERATURE (JOULES) 
PLASMA TEMPERATURE (JOULES) 
ARC CURRENT 


94 


CC? CC Oa) a) C2 Ce Grae 


ND 

DTK 
COULOG 
ECODE 


Te 


KK KR KK KR Kee ee 


* 


NUMBER OF DATA SETS 
CATHODE TEMPERATURE INCREMENT 
COULOMB LOGARITHM 


ERROR CODE 
O= GOOD 
l= ERROR 


ieoeb Or BMS nN 
l= THERMIONIC 
2= INTERMEDIATE 
3= FIELD 
4= NONE 


SUBROUTINES REQUIRED 
BOMCh PHO, .CODE, TYEa) 


SKORIK OOK OK OIOK OK KOK OOK OK OK OK IK OOK OOK OK KOK OK OK KOK 


ies Oye Ow, 2 1,P HI BRT A MARR, Z,NO,M,E,K,H,EPS 
RE oun Di RA! arr lene olG ena Vi VRIVK2, VT 
REO ol, 02,05,04,05,0,Jn,J1, KIRKE Din cOuLOG 
INTEGER ECODE,TYP,N,ND 

OPEN(6,FILE="ARC.DAT ,STATUS="NEW’) 


OPEN(7, FILE="PRN’) 


WRITE(*,*) “ENTER 
READ(*,*) ND 
WRITE(*,*) “ENTER 
READ(*,*) KTP 
WRITE(*,*) “ENTER 
READ(*,*) TK 
WRITE(*,*) “ENTER 
READ(*,*) DTK 
WRITE(*,*) “ENTER 
READ(*,*) I 
Whelieie, 000) ND, 
WRITE(7,900) ND,I 
MI=63.54 
LO=300300. 
Bi=2o//D-18 
PHIO=4.4 

AI=40. 

BETA=3.21 

Z=1.8 
NO=6.0220D23 

B= saltso D3 1 
B=LoOZ2-19 
K=13807D-23 
HeibroZzocp~34 
PI=3.141592654 
EPS=8.8542D-12 
BER=L0 /NO 

l= 0 DA 


NUMBER OF DATA SETS: ~ 

INITIAL PLASMA TEMPERATURE (EV): ° 
INITIAL CATHODE TEMPERATURE (K): ° 
CATHODE TEMPERATURE INCREMENT (K): ~ 


ARC CURRENT (A): ° 


SIs 


S2ZSa07 1) 
S3=0" 
Saal gael 
55=153.32236e420 
Vig=ibe 
MI=MI*1L6726D-27 
KT P= Xe 
DO 20 N=2ND 
KTK=K*xTK 
JES Or 
PTK=(10**(S1/TK +S52*DLOG10(TK)+S3*TK+S4))*S5 
10 ATP=1/DSQRT(1l. + PTK/KTP*(2*PIKH**2 /(M*XKTP))**(3./2.)* 
1 DEXP(EI/KT P)) 
ETA=(8./3. +ATP-DSQRT(5.*PI*XKTK/(6.*KTP)))/ 
di (8./3. +ATP+DSORT(S.4 FIG. 
JIZE*PTK*XETAXATP/DSQRT(2* PIXMI*XKTK) 
F=DSQRT(JI*DSQRT(8* VK*MI/E)/EPS *(1-JE /JIXDSQRT(M/MD)) 
PHI= PHIO*E-DSQRT(E**3*EF /(4*PIXEPS)) 
CALL EMM(TK,F,PHIO, JE, BCODESTYP) 
R=DSQRT(I/PI*(JE+JD)) 
VKi=((JI/JE)*(EI+2*KTP)+BETA*KTPxX(1.4+JI/JE)-2*KTK)/E 
G=PTK*(1-ETA)/DSQRT(2* PIXMI*XKTK) 
VK2=4*XIXTK/(PIXR¥*JDI+(JE/JI)*(PHI+2 *KTK)+(G/JI)XEER*E 
1 -EI-2*KTP+PHI)/E 
IF(DABS(VK2-VK1i).GT. LOD-1) THEN 
IF(VK2.GT.VK1i) THEN 
VK=VK1i 
IF(DABS(VK2-VK1).GT. 100) THEN 
KTP=KTP+iQGeak 
ELSE 
KTP=KTP+DABS(VK2- VK1)*K 
ENDIF 
ELSE 
VK=VK2 
IF(DABS(VK2-VK1i).GT. 100) THEN 
KTP=KTP-101*K 
ELSE 
KTP=KTP- DABS(VK2- VK1)*K 
ENDIF 
ENDIF 
IF(VK.LT.0) VK=DABS(VK) 
GOTO 10 
ENDIF 
VK=(VK1+VK2)/2 
J=JE+JI 
COULOG=23.0-DLOG(DSQRT(PTK*ATPXETA/K)/AKTP/K)**2) 
SIG=L53D-2*(KTP/K)*x«(3./2.)/COULOG 
U=(IKE*GK(2* PI*SIG*KR*J))*K(L+Z)/L +(E*XG/J)*(L +Z))- 
1 BETAX*(KTP/E)/A1L+(E*G/J)*(1.+2Z)) 
Vi= Via 
WRITE(6,9 10) TKSKTE Ey oceys 


96 


900 


WRITE(6,920) PTK/1.01D5,ATP,ETA,JI/J,R 
WRITE(6,930) F,JE,JLJ 

WRITE(6,940) COULOG, SIG, G 

WRITE(7,910) TK,KTP/E, VK,U,VT 

WRITE(7,920) PTK/1.01D5,ATP,ETA.JI/J,R 
WRITE(7,930) F,JE,JLJ 

WRITE(7,940) COULOG,SIG,G 

Gr=TK+DTK 

CONTINUE 

STOP 

FORMAT(1X,13,5X,F 4.0) 
FORMAT(1X,F5.0,5X,F6.2,5X,F6.2,5X,F6.2,5X,F6.2) 
FORMAT(1X,F4.1,5X%,F6.4,5X,F6.4,5X,F6.4,5X,1P, D10.4) 
FORMAT(1X, 1P, D10.4,5X,D10.4,5X,D10.4,5X,D10.4) 
FORMAT(1X,F4.1,5X,F7.1,5X,1P, D10.4) 

END 


mei 
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